Field theoretical study of a spin- 1/2 ladder with unequal chain exchanges 
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We study the low-energy properties of a Heisenberg spin- 
1/2 zigzag ladder with different exchange constants on the two 
chains. Using a nonlinear cr-model field theory and abelian 
bosonization, we find that the excitations are gapless, with a 
finite spin wave velocity, if the values of the chain exchanges 
are small. If the chain exchanges are large, the system is 
gapped, and the energy spectra of the kink and antikink ex- 
citations are different from each other. 

PACS numbers: 75.10.Jm, 75.10.Ee, 71.10.Hf 



I. INTRODUCTION 

For the last several decades, one-dimensional and 
quasi-one-dimensional quantum spin systems have been 
studied extensively due to their unusual properties. Ex- 
perimentally, many such systems are known to present a 
wide range of unusual properties, and a variety of analyti- 
cal and numerical techniques exist for studying these sys- 
tems theoretically. The three observations which make 
the low-dimensional spin systems particularly intersting 
are, (i) Haldane's conjecture for one-dimensional antifer- 
romagnetic spin systems , (ii) the discovery of high- 
temperature superconductivity and its magnetic proper- 
ties at low doping ||, and (iii) the discovery of ladder 
materials 

In spin ladders, two or more one-dimensional spin 
chains interact with each other. For ladders with the 
railroad geometry, it has been observed that spin-1/2 
systems with an even number of legs are gapped, while 
systems with an odd number of chains have gapless ex- 
citations However, the frustrated zigzag ladder || 
shows gapless spin liquid state or the gapped dimer state, 
depending on the ratio of the exchanges of the rungs to 
the chains pQ |. 

Another interesting kind of system is the spin-Peierls 
system such as CuGeO^. The authors of Ref. (s) have 
explained the form of the ground state and the low- 
temperature thermodynamic properties of this sample by 
showing that there is spontaneous dimerization of the 
nearest neighbor interaction below a particular tempera- 
ture. Their model includes a dimerization in the nearest 
neighbor exchange coupling. 

Relatively less is known about a spin ladder with asym- 
metry in the chains. An extreme case of this situation, 
i.e., when a chain is absent from the Fig. 1 (sawtooth 
chain), has been studied by two different groups 0,0. 
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FIG. 1. Schematic diagram of a spin ladder with unequal 
chain exchanges. 



The ground state is like that of the Majumdar- Ghosh 
model |l3| , except that the two kinds of low-energy exci- 
tations (kinks and antikinks) which interpolate between 
the two degenerate ground states have quite different ex- 
citation spectra. 

In this paper, we study a two-chain ladder system with 
unequal exchange constants of the chains, rather than a 
dimerization in the rungs. Very recently, Voit et al H] 
have studied this problem. We will discuss their results 
below. The plan of the paper is as follows. In Sec. II, 
we will analyze the problem using the nonlinear cr-model 
(NLSM) field theory. Sec. Ill will discuss the abelian 
bosonization approach. 



II. NONLINEAR a-MODEL STUDY 

In this section we study the NLSM field theory of our 
model. The schematic plot is shown in Fig. 1. The 
system can be viewed either as a ladder with unequal 
exchanges on the two chains, or as a single chain with 
unequal next-nearest neighbor exchanges. The Hamilto- 
nian is 



H 



E 



[JA-5 n+1 + J 2 (l-(-l) n S)S n -S, 



n+2 



(1) 



where n is the site index, Ji, J2 > 0, and < 8 < 1. If 
we view the system as two chains, then the lower (upper) 
chain contains the odd (even) numbered sites. 

Using the NLSM field theory, one can describe the 
low-energy and long-wavelength excitations. In the Neel 
phase, this is given by an 0(3) NLSM with a topological 
term. Here we define two fields <f> n and l n as a linear 
combination of two spins as 
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I, 



<J2n-l ~~ S211 
2S 

02n-l + S211 

2a 



(2) 



where a is the lattice spacing. It can be easily checked 
that 



• L = 



1 



1 

S 



2Jtl 
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Thus 4> n becomes a unit vector in the large S limit. The 
unit cell of the classical ground state of the Neel phase is 
labeled by an integer n, and it contains the sites 2n — 1 
and 2n respectively; the length of a unit cell is 2a. 
The fields 4> n and l n satisfy the commutation relations 
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where m and n are the unit cell labels, a, b, c denote the 
x, y, z components of the field, and e xyz is the completely 
antisymmetric tensor with e xyz = 1. This relation en- 
ables us to write l n = 4> n x n n , where the vector II is 
canonically conjugate to <p namely, 



2a 



b • 



(5) 



To define the continuum limit of this theory, we introduce 
a spatial coordinate x which is equal to 2na at the loca- 
tion of the n th unit cell. Summations are then replaced 
by integrals, i.e., ~ ¥ I dx/(2a). 

Since <f> n is a unit vector, (f> and <fi' are orthogonal to <f>. 
In the low-energy and long- wavelength limit, the dom- 
inant terms in the Hamiltonian are those which have 
second order space-time derivatives of <j) and first order 
derivatives of I (since I contains the the first order deriva- 
tives of 4>) p5t^6| . To get the continuum Hamiltonian, we 
expand the fields 4> n +\ = 4>{x + 2a) and l n +\ = l(x + 2a), 
where x = 2na, as 

$(x + 2a) = (f{x) + 2a$(x) + 2a 2 $"(x) + ■ ■ ■ , 
l(x + 2a) = l(x) + 2al'{x) + ■■■ . (6) 

Using Eqs. (||), (0), and @, we obtain the Hamiltonian 

- 2 fi - . r 

, (?) 



H = ] d x { C JL {l+ _J ? + _ 



where 



= 2J 1 Sa^l - U2/J1 , 
2 



Sy/1 - U2IJ1 ' 

= 2tt5 . 



Note that the values of c, g 2 and 9 turn out to be inde- 
pendent of 8 in this approach. Further, this NLSM is 
valid only if J 2 /Ji < 1/4. For J 2 /Ji > 1/4, a different 
NLSM is required. 

One can find the energy-momentum dispersion relation 
of the form u) — c\k\, where c is the spin wave velocity, 
by considering small fluctuations around <j> — (0,0,1), 
and expanding the Hamiltonian in (0) to second order in 
those fluctuations. Similarly one can find the strength of 
the interaction between the spin waves, g 2 , by expanding 
the Hamiltonian to fourth order in the fluctuations. 

From Eq. (||), we see that the coefficient of the topo- 
logical term 9 — tt for a spin-1/2 system. This implies 
that there is no gap in the low-energy excitation spec- 
trum. This result is different from that of the zigzag 
ladder with a dimerization in the nearest neighbor inter- 
action, i.e., with a term like J\{— l) n S n ■ SVi+i- In that 
case 9 — 2ttS (1 — ft) is different from tt for a spin-1/2 
system |l6|,[l7]], and the the low-energy excitations are 
gapped. 



III. ABELIAN BOSONIZATION STUDY 

In this section we study the low-energy excitations 
spectrum using abelian bosonization. We express our 
Hamiltonian as the following sum, 



where 



H = Hi + H 2 + H 2S , 
Hi = Ji S n ■ Sn+i j 

a 

Hi = J2 J]] S n ■ S n +2 , 

n 

H 2 8 =-J 2 Sj2 (- 1 )" Sn ■ S n+2 . 



(9) 



(10) 



(8) 



We can convert this Hamiltonian to a Hamiltonian of 
spinless fcrmions through the Jordan- Wigner transfor- 
mation. The relations between the spin and the electron 
creation and annihilation operators are, 

S* = fy n - 1/2 , 

n-l 

S n = (~l)"^n exp[i7T ^ n j\ > 

i=— 00 

n-l 

i=— 00 

where rij — ipjipj is the fermion number at site j. The 
Hamiltonians in (|l0|) then become 

Hi =-y Wt+l'&n + *l>bl>n+l) 

n 

+Ji - l/2)(^ +1 Vn+i - 1/2) , (12) 
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H 2 =J 2 J2 Wl+aV-n + h.c.)(^+i^n+i - 1/2) 

n 

+ J 2 J2 - i/mi +2 ^ n+2 - 1/2) , (13) 

■n 

H 2S = -J 2 ^(-l)"(^ +2 Vn + ll.C.)(^ +1 ^n+l - 1/2) 
n 

-J 2§ J2(-m^n - 1/2)^+2^+2 - 1/2). 

n 

(14) 

We will assume below that J\ >> J 2 . Since there is 
no applied magnetic field, the two Fermi points occur 
at kp — ±7r/2. We can linearize the energy spectrum 
around these Fermi points, and express the lattice oper- 
ators in terms of two continuum fields R and L which 
vary slowly on the scale of a lattice spacing, 

1> n = [ i n R(n) + (-i) n L(n) ] , (15) 

where R and L are describe the second-quantized fields 
of right- and left-moving fermions respectively. Now we 
bosonize our Hamiltonian following the standard proce- 
dure. The basic relations used to obtained the bosonized 
Hamiltonian are as follows pq . 

S*(x) = a [ p(x) + (-iy M(x) ] , (16) 

where the fermion density p(x) —: R^ (x)R(x) : + : 
L'(x)L(x) :, and the mass operator M{x) —: R^{x)L(x) : 
+ : L^{x)R(x) :. (The double dots denote normal order- 
ing). The bosonized expressions for p and M are given 
by 

p(x) = - —= d x (f>(x) , 

M{x) = — cos(2^(a;)) . (17) 
The bosonized version of H\ is known to be 
Hl =jdx[ V -±ftf + ^ M f 

+ : cos(4^) : ] , (18) 

where v\ = irJ\a/2 is the spin wave velocity, and K 
1 1 is the bosonization interaction parameter for the 
isotropic spin- 1/2 antiferromagnet. The last term in Jl8| ) 
is marginal, and is known to have no effect at long dis- 
tances in the sense of the renormalization group (RG) 
fjj|. Using Eqs. (0), (|l|), and ©, we get the follow- 
ing expression for H 2 , where we have ignored terms of 
the order of a 4 and higher, 

H 2 = J 2 a 2 h (Pn+i ~ (~l)"A/„ +1 ) x 

n 

(Pn + Pn+ 2 + (-1)"M„ + (-l) n M n+2 ) 

+ (p n + {-l) n M n ){p n+2 + (-l) n M n+2 )}. (19) 



To derive this expression, we have used Taylor expres- 
sions such as 

R(n + 2) = R(n) + 2aR'(n) + 2a 2 R"(n) + ■•■ 

(20) 

to write 

R\n + 2)R{n) + R\n)R{n + 2) 
= R\n + 2)R{n + 2) + R\n)R{n) + 0(a 2 ) . (21) 

On keeping only the terms which do not oscillate as (— l) n 
(which would give zero in the continuum limit) and then 
expressing the operators p and M in the bosonic lan- 
guage, the above expression becomes 

H 2 =v 2 [ dx [--{d x (j>) 2 + — ^jB»(4^)], (22) 
J ir 2(7raj z 

where v 2 — J 2 a. Both the terms in ( |22] ) have scaling 
dimension 2, and they are marginal. It is known that 
they have no effect in the RG sense as long as J 2 /Ji < 
0.241 |ist| . Finally, one can find the bosonized expression 
for H 2 $. To begin with, the nonoscillatory part of the 
Hamiltonian is given by 

H 2 s = -J 2 Sa 2 ^2 [p n M n+2 + M n p n+2 - 

n 

Pn+1 (M n + M n+2 ) + M n+1 ( Pn + p n+2 )}. (23) 

Now we perform an operator product expansion of the 
above Hamiltonian. In the limit z — > w, we can use the 
following expansion ]20| , pT[ |, 

dMz) : e 1 ^ : = ^- : e ip ^ : 

z — w 

+ : d z <i>{z)e lfi ^ z) : (24) 

for K = 1/2. The second term in ( pi| ) is a total deriva- 
tive, and its contribution will therefore vanish in the 
Hamiltonian where it appears inside an integral over all 
x. From the first term in Eq. (^), we see that the var- 
ious terms in Eq. (^3|) cancel each other in pairs which 
have z — w — ±a or ±2a. If this cancellation had not 
occurred, the bosonized version of H 2 $ would have been 
proportional to the operator cos(2y / 7t0), which has scal- 
ing dimension 1 and would therefore have been relevant. 
However, due to the cancellation, H 2 s contains no rele- 
vant operators with scaling dimension less than 2. Thus 
the system continues to remain gapless (and lies in the 
same spin-liquid phase as the model described by the 
nearest neighbor Hamiltonian Hi) even if 6 ^ 0, provided 
that J 2 << J\. This is in contrast to a dimerization in 
J\\ in that case abelian bosonization correctly produces 
a relevant term which leads to a gapped phase. 

Recently Voit et al |l4| have studied the same model 
as ours using abelian bosonization followed by a pertur- 
bative renormalization group analysis. They claim that 
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the perturbation due to the unequal chain exchanges ( J 2 <5 
terms) is relevant but that it does not lead to a gapped 
phase; instead, they argue that it leads to a different 
fixed point where the system is gapless and has a vanish- 
ing spin velocity. In our study, we have shown that due 
to some cancellations, the J 2 <5 terms do not lead to any 
relevant terms in the continuum theory, and therefore 
that the system remains gapless if J 2 is small. 

Finally, let us discuss the low-energy excitations of this 
model when J 2 becomes larger. In particular, we find 
that these are given by kinks and antikinks when J 2 = 
Ji/2, for all values of 5 lying in the range [0,1]. (The 
Majumdar-Ghosh model is a special case of this where 
(5 = 0). Then the Hamiltonian in (|l[) can be written, up 
to a constant, as the sum 

H = J2 [( 1 + ^)(S ? 2„-i+5 2n + 5 2 „ +1 ) 2 

n 

+ (1 - S)(S 2n + S 2n +i + SW2) 2 ] • (25) 

Hence the ground state is given by a configuration in 
which the total spin of each triangle is 1/2. Since this can 
be done either by forming a singlet with the pair of spins 
(2n — 1, 2n) for all values of n, or by forming a singlet with 
the pair of spins (2n, 2n + 1) for all values of n, we see 
that the ground state is doubly degenerate. Let us denote 
these two ground states by A and B respectively. The 
lowest-energy excitations (kink and antikink) are formed 
by interpolating between these two ground states. The 
kink has the ground state A on the left and the ground 
state B on the right, while the antikink has the ground 
state B on the left and the ground state A on the right. 
For general values of 5, we find that the kink and an- 
tikink dispersion are nondegenerate in contrast with the 
Majumdar-Ghosh model ]ll] , p^| . A simple variational 
calculations gives the kink and antikink dispersions to 
be Ji(l - <5)(5 + 4cosfc)/8 and Ji(l + <5)(5 + 4cosfc)/8 
respectively. Hence the minimum gap for the kink and 
antikink excitations are Ji(l — S)/8 and Ji(l + <5)/8 re- 
spectively, occurring at k = it. This result has been also 
obtained by the Voit et al Jl4| . 

To summarize, we have studied the low-lying excita- 
tions of a zigzag ladder with unequal chain exchanges. 
Both the NLSM and abelian bosonization show that the 
system remains gapless if J2 < < J\ . We have also shown 
that the system is gapped (with two degenerate ground 
states) for J 2 = Ji/2. It would be interesting to use 
numerical techniques like the density-matrix renormal- 
ization group ^2j, to study the complete phase diagram 
of the ground state as a function of the two parameters 
J 2 / J\ and S. 

We (SS and DS) would like to thank Marco Ameduri 
for useful discussions. DS would like to thank Bruce Nor- 
mand for helpful comments. 
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